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Abstract 

Extensions of the Standard Model are being considered as viable settings 
for a first-order electroweak phase transition which satisfy Sakharov's three 
conditions for the generation of the baryon asymmetry of the Universe. These 
extensions provide a sufficiently strong phase transition and remove the main 
obstacles which appear in the context of the Standard Model: A far-too- 
high lower bound on the Higgs mass, immediate wipeout of the newly-created 
baryon asymmetry, and insufficient CP violation. We describe the Universe 
hydrodynamically as a fluid coupled to the Higgs field via a phenomenological 
friction term, and study the time evolution of bubbles nucleated during the 
phase transition. We express the friction term in the hydrodynamic equations 
in terms of the particle content of the model, calibrate the friction on the basis 
of existing calculations for the Standard Model, and produce predictions for 
the velocity of the expanding bubble wall in the stationary regime. This 
way we develop a very efficient approach to compute bubble velocities. As 
an example, we apply our formalism to the first-order phase transition of 
a dimension-6 extension of the Standard Model which, within the present 
bounds on the Higgs mass, can reproduce the observed baryon asymmetry 
of the Universe. Depending on the strength of the phase transition, the wall 
velocity varies from about 0.3 to approaching the speed of light. Our method 
can easily be adapted to compute wall velocities in other interesting extensions 
of the Standard Model. 
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1 Introduction 



A first-order electroweak phase transition (EWPT) in which two differentiated phases 
briefiy coexisted could be the setting for the genesis of the observed baryon asymme- 
try of the Universe [Il[2l[3lll] (for a review, see e.g. [5]), and may have left signatures, 
like a gravitational wave background or primordial magnetic fields, which we might 
be able to observe in relatively short order [SI El El El HOI [H]- Such a transition may 
satisfy Andrei Sakharov's three conditions for electroweak baryogenesis QJ: CP vio- 
lation, baryon number violation, and deviation from thermal equilibrium, the latter 
imposed by the expanding bubbles of the broken symmetry phase on the whole of 
the medium at one point or another until the time when the whole Universe has 
crossed over to the new phase and the transition ends. A first-order EWPT requires 
the presence of two stable or metastable phases which correspond to extrema in the 
high-temperature expansion of the Higgs thermal effective potential (Figure [1]). It 
is well known that the observed baryon asymmetry could not have been produced 
in such a scenario within the Standard Model (SM) unless the Higgs mass were far 
lower than allowed by present experimental boundqj. This is because sphaleron tran- 
sitions must be sufficiently suppressed inside the bubbles of the new phase to avoid 
washout of the newly-created baryon asymmetry [T6l fT7]. Also, the electroweak 
phase transition is not even first-order in the Standard Model except (again) for 
Higgs masses far below present experimental bounds [IB]. Simple extensions of the 
standard model, however, are still capable of producing a sufficiently high baryon 
asymmetry through a first-order EWPT for presently acceptable values of the Higgs 
mass m//, be it via extra Higgs fields [ISl EHl EI] , light top squarks [221I231IM 125] . 
singlets [ll[26l|27l|28l|29l|30l|3l]or non-standard Higgs potentials [32l [33]. Below 
we will consider the latter possibility. 

Calculations of the dynamical characteristics of a first-order phase transition 
(like the expansion velocity of the walls of bubbles of the new phase after nucle- 
ation) mostly rely on solving the relevant hydrodynamical equations, modelling the 
radiation-dominated early universe as a perfect relativistic fluid [311 ESI EHl EZl EH 
ES So] . Usually, the WKB (semiclassical) approximation is adopted. This assumes 
that the de Broglie wavelength of the particles in the plasma is considerably smaller 
than the width of the bubble wall, which should be true for all particles in most 
situations except possibly for very infrared bosons, which are assumed to have little 
influence on the wall dynamics [36] • Such treatments derive a fluid equation which 
describes the dynamics of the variables characterising the cosmic fluid (fluid veloc- 
ity, temperature, and Higgs vacuum expectation value). The fluid equation includes 
a friction term which tends to slow down the propagation of the bubble wall. In 
some treatments a phenomenological parameter which characterises the medium is 
introduced into such a term [351 HQ]. In others an explicit calculation of the friction 

■^The Atlas and CMS experiments at the Large Hadron Collider (LHC) had bracketed the Higgs 
mass to between 117 and 127 GeV when the discovery of a 'previously unknown boson', consistent 
with the Higgs and with a mass of 125 — 127 GeV, was announced in July 2012 [T2l [T3l fT4] 
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from the modelling of the plasma is carried out [511 EH]- Friction is a function of 
the deviation from thermal equilibrium of the particle components of the plasma 
which couple strongly to the background Higgs field. If local thermal equilibrium 
(LTE) were maintained throughout the plasma there would be no friction and the 
expanding bubble walls would accelerate to the speed of light [H]. The presence 
of friction may lead to two outcomes: Either the advancing bubble wall reaches a 
steady state under the influence of friction and propagates with constant velocity, 
or friction is too weak and the wall accelerates without bound towards ultrarela- 
tivistic velocities. The feasibility of such 'runaway walls' has recently been proven 
theoretically [32|[39]. 

In the following we lay out a general approach to bubble velocities, based on an 
efficient parametrization of friction exerted by the plasma on the bubble wall. In 
ref. |12], we have used this framework for the first time to compute the wall velocity 
in the MSSM for realistic Higgs masses, requiring use of the 2-loop thermal Higgs 
potential. Here we take the example of the dimension-6 extension of the Standard 
Model with a cut-off M, discussed in detail in refs. |33], to demonstrate the power 
of our approach to bubble velocities. This extension if the SM provides also new 
sources of CP-violation necessary to obtain realistic values of the baryon asymmetry 
within the model parameters |43j and results in a phase transition strong enough 
to avoid washout. We apply the semiclassical, perfect-fluid approximation to the 
electroweak plasma with our extended thermal potential. 

We adapt standard hydrodynamical calculations to our model which employ a 
phenomenological friction parameter. In order to calibrate this parameter realisti- 
cally for our model we base ourselves on existing results for the wall velocity within 
the Standard Model and the description of the evolution of particle populations 
through Boltzmann evolution equations in the form known as the relaxation time 
approximation (SB]- We compare the calibration from Standard Model results with 
that found through the theoretical criterion for 'runaway walls' mentioned above. 
Thus we predict values of the friction parameter with which the steady-state ex- 
pansion velocity of the bubble walls may be calculated for a range of the model 
parameters M, niH- 

We flnd the wall velocity to be often subsonic (lower than the speed of sound 
in the medium), similar to existing results for the minimal SM. For relatively low 
values of the Higgs mass and the dimension-6 cutoff scale M, however, the phase 
transition becomes very strong and the steady-state wall velocity may shift into 
the supersonic regime. Eventually we predict runaway, rather than steady-state, 
behaviour for the highest values of the strength of the phase transition within our 
parameter space. We comment briefly on the possibility of 'hybrid' steady-state 
solutions in the frontier between the subsonic and supersonic regimes, as well as on 
the existing criteria for stabihty of the growing bubbles as they apply to our model. 

Lastly, we provide a 'user-friendly' description of our friction parameter cali- 
bration, based on the relaxation time approximation, for use in any model given a 
reference calibration point for the friction. The results presented here can directly 
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be used to treat cases with SM-like friction, e.g. SUSY scenarios without hght stops, 
as for instance can be reahsed in singlet extensions of the MSSM. 

2 Single field hydrodynamics and the electroweak 
phase transition 

2.1 The effective potential and critical free energy 

The fundamental tool employed to study the electroweak phase transition is the 
finite-temperature effective potential for the background Higgs VEV 0, V{(j),T), 
representing the free energy per unit volume of the field configuration (see e.g. 
At high temperature, the thermal potential has only one local minimum (Figure [1]) 
at = 0, representing the unbroken symmetry phase of the theory and the state of 
the universe before the phase transition. As temperature decreases a second local 
minimum of the potential appears for a 0o > 0, corresponding to the broken sym- 
metry phase of the theory. As at first ^(00)^) > ^(0,T) the broken symmetry 
phase is energetically disfavoured and the phase transition may not yet proceed. As 
temperature decreases further, V{(j)o,T) also decreases until it becomes degenerate 
with the 0-VEV minimum at the critical temperature T^. From that point on con- 
figurations ('bubbles') of the new phase may, in principle, form and the transition 
may take place. 

The free energy of a bubble of the broken symmetry phase is given by 

F = l^d'x[^{V<pr + V{<P)]. (1) 

The first term represents the positive contribution from the bubble wall, defined 
as the region where the Higgs VEV varies between its zero value in the symmetric 
phase outside the bubble and its nonzero value inside. The second term represents 
the negative contribution from the value of the Higgs potential in the interior volume 
of the bubble. Assuming spherically symmetric bubbles, for each temperature there 
exists a critical radius such that nucleated bubbles larger than the critical size 
will expand spontaneously (thus minimising their free energy) whereas smaller than 
critical bubbles will collapse under surface tension. The critical radius is infinite at 
Tc and decreases monotonically with temperature. The critical bubble is a static 
solution of the Euler-Lagrange equation 



^ 6L 6L 

In spherical coordinates with radial coordinate r (and since for a static solution 
dt = 0) the Euler-Lagrange equation becomes 
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dr'^ r dr 

which has to be solved with boundary conditions 



(j){+oo) = 
= 



+ zl = l^ (3) 



dr 



J r=0 



The critical bubble solution can be integrated numerically (eq ([T])) over the 
bubble volume to find the critical free energy. The radius of the critical bubble at 
the critical temperature becomes infinite and the equation for the Higgs VEV may 
be written as 



d^^dV 

dz^ d(f) ^ ' 



with boundary conditions 



0(-oo) = 

0( + Oo) = Vc 

where Vc is the nonzero Higgs VEV at the degenerate minimum at T = (Fig [2]). 

2.2 Nucleation and finalisation temperatures 

The probability of bubble nucleation within the horizon volume in the time interval 
dt may be calculated from the critical free energy as (see again e.g. [1]) 

dP={TlV).Vn-dt=^^e-^^l^^^ 

where T/V= A^(T)e~^'=/^ ~ 2^4g-Fc/r jg ^^^q nucleation rate per unit volume with 
Fc the critical free energy at the temperature T. The horizon volume is Vh = d^^j = 
i2tf = H-'^ with H-^ = and Mpi the Planck mass. ^ = (1/34) at the time 

of the electroweak phase transition. The phase transition starts at the nucleation 
temperature for which 

P{T = T„) = / dP = l (5) 
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A true-vacuum (broken phase) bubble nucleated at a temperature T' and ex- 
panding with constant velocity (3 will have grown to a radiu^ at a later temperature 
T 



(6) 



Thus the fraction / of the causal volume which has crossed over to the broken 
phase at a temperature T < T„ is given by (neglecting mergers between bubbles) 

1 in 



1 Att = 

^(^) = / RB{TXfdP{T') 



3 Jt \ ^ T J \T T'J H{T'Y V 



3 ' Jt V T' J T'5 

f{Tf) = 1 defines the phase transition's finalisation temperature at which the tran- 
sition can be assumed to end. 

2.3 Equation of motion for the Higgs VEV 

The equation of motion for the background Higgs field (see e.g. ^J) is a function 
of the phase space population density functions f{p,x) for all the particles present 
in the plasma, 

where V{(f)) is the renormalised vacuum potential and the sum is over massive degrees 
of freedom so that, for the appropriate couplings, e.g. = for fermions, 
_ bosons. The population densities can be expressed as a thermal 

equilibrium part plus a deviation, / = /q + 6f, with the equilibrium distribution 
for fermions/bosons, in the rest frame of the plasma, given by foiT) = ^e/t^i (with 

E = \J + IpI^). The derivative of the vacuum potential combines with the integral 
of the equilibrium part of the distributions to give the finite-temperature effective 
potential V^(0, T). The equation of motion then becomes 



^We ignore the effect of the expansion of the Universe, since the transition is extremely fast, so 
that dR = (3dt. 
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The (5/-dependent term is the friction term. Thus, as stated earher, friction effects 
are the result of the deviation of the particle populations in the plasma from equi- 
librium. If these did not exist there would be no friction (and the bubble wall would 
propagate at the speed of light). 



2.4 The pressure on the wall 

A direct application of eq (|9]) is the calculation of the pressure felt by the wall |36] . 
Assuming a planar wall advancing with a steady velocity in the z direction we may 
rewrite (jH]) in the rest frame of the advancing wall as 

Multiplying the whole expression by ^ = 0' and integrating over z across the wall, 
we find the pressure on the wall. Note that at both ends of the wall ^ = 0. 

Therefore I (h"(h'dz = =0. We use the fact that the finite-temperature 

L ^ Jo 

effective potential of the Higgs field equals the free energy density per unit volume 
F, which is equal to minus the pressure. Thus integrating along the 2;-direction 
yields 

In a steady state, the pressure felt by the advancing wall equalises with the integral 
of the friction term as written in eq. flTT]) . 



2.5 Boltzmann evolution equations 

In order to solve the equation of motion for the Higgs field, 6f must be expressed 
in an appropriate form. Within the WKB approximation the evolution of the par- 
ticle populations is usually accounted for through the relevant Boltzmann evolution 
equations, 

I = 9,/ + f ■ d,f + fl- d^yf = -C[f] (12) 
where C[f] is the so-called collision integral, written in its full form as [36] 



xi27rf5\p + k-p'-k')Pm (13) 

and P[/J = /i/2(l ± /3)(1 ± h) - hUil ± /i)(l ± /2)- The sum is over all relevant 
four leg scattering diagrams, p and k are the incoming, p' and k' the outgoing 
momenta. M is the scattering amplitude for each process. The fi are population 



8 



factors. The first (positive) contribution to P represents a particle with momentum p 
being removed from the state by a colhsion, weighted by the populations of colliding 
particles. The second (negative) contribution represents a particle being scattered 
into the state. The (1 ± /) factors stem from particle statistics (— for fermions and 
+ for bosons). 

The equilibrium distribution for fermions/bosons can be generalised to a frame 
in which the fluid is not at rest but has a bulk velocity v in, say, the z direction. 
The equilibrium distribution becomes 



/o 



1 



(14) 



with /3 = ^ and 7 = :^7Y=f • Taking v = recovers the version for a fluid at rest. 
In order to work with the Boltzmann equations we need the relevant position and 
momentum derivatives of the equilibrium distributions. 



dfo 
dpz 



(^(.M{E-vp.) ± 1)2 



(15) 



and 



dz 



(16) 



The hydrodynamical calculations across the bubble wall profile are usually car- 
ried out in the rest frame of the (steady-state) bubble wall, in which dt = for all 
quantities. A planar wall will 'see' the plasma as moving with a given 'fluid' velocity. 
Note that, if the plasma itself has a position-dependent, nonzero bulk motion, in- 
duced by the passing of the wall, the plasma velocity 'seen' by the wall at any given 
position may not coincide with the steady-state velocity of wall propagation. Taking 
dp6f = 0, and noting that Vz ■ dz6f = ^ ■ dz6f and pi = —^Uz ■ 



1 d(m^) 



Boltzmann equation in the rest frame of the advancing wall becomes 



the 



{m^)' pz 
2E 



- V 



Pz [m 



E 2E 



[m 



2V 



2E 



-vf3'y 



dtf + X ■ d^f + p- djjf 



Pz 

E 



dzifo + 5f) 



[m 



2V 



2E 



-dp Jo 



2\n ^l3-i{E-vp,) 

(^^P^{E-vp,) ±1)2 + 'E^'^^f = 

'^M(E-vp.) 1)2 + '^^f = -'^[f]- 



(17) 



2.6 The relaxation time approximation 

A number of simplifications are possible for the collision integral C[f]. The free 
particle approximation jUJ US] assumes C[f] = 0. This may represent the case in 
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which particle free paths are much larger than the thickness of the wall. In this 
approximation the Boltzmann equations can be solved exactly (see e.g. [SS])- 

Our starting point will be the relaxation time approximation in which we assume 
C[f] = where the relaxation timescale r is usually considered independent of 
momentum [H]. Following this approximation the Boltzmann equation becomes 

-^/577-wi^Z7^r^^ + T^'^/ = — • (18) 



2E ' {eP-<iE-vp.) ±lf E ' T 

The usual procedure is to neglect either the 6f or the 6f' term. The authors of 
e.g. |36] assume L ^ t for the wall thickness L, so that Sf ^ ^ This is a 

natural assumption if the bubble wall is expected to be relatively slow since slower 
walls are thicker and not as sharp (quantities vary more slowly across the wall) and 
relativistic time dilation is not a factor in increasing the characteristic timescale 
for particle interaction (the relaxation time). In this approximation we drop the 
derivative term in f|T7|) and obtain (for one degree of freedom), 

(^2y g/37(£-!;p^) 

Inserting this form for 6f into the friction term in the equation of motion for the 
background Higgs field ( ITU]) we obtain 

~d^ J (2^0^ ~ 
dm^ r d^p {my e^T(^-^P^) 

~d^ J (27r)32E^^^^ 2E (e/37(i?-«P.) ± 1)2 (^^^ 

which gives an equation of motion for a slow wall, in the rest frame of the advancing 
wall, of the form 

dV(6,T) ^dm^ f d^p ^ (m^ e^i{E-vp.) 



dz^ d4> ^ d4> J (27r)32E ' 2E {ef^^(E-vp.) ± ly 

The pressure difference felt by the advancing wall (for one degree of freedom) in this 
approximation is 

AP3W = j dz^ j j^^r^iv 2^ (g,,(^_.,,)^,),. (22) 

The opposite limit is that of a fast, relativistic wall. This will naturally be a thin, 
sharp wall and relativistic time dilation will lenghten the characteristic relaxation 
time for particle populations. Thus we assume L <^t and 5f' ^ ^ ^ Therefore 
in this limit we may neglect the 5f term in (|T8|) and obtain 

p (rn^y qPi{e-vp-^) 
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which (assuming for simphcity E ^ in the exponentials) can be integrated to 



^2 ^^^{E-vp,) 



The friction term becomes in this hmit 

dm^ f d^p 



{2'Kf2E 



dm^ f d^p ^ e^i{E-vp.) 

-vp'j- 



and the equation of motion for the fast wall case 



(25) 



ciV dV{(l),T) ^dm^ f d^p e^i{E-vp.) 



dz^ d(f) ^ d(f) J {27T f2E2p, ^ ' [e^iiE-vp.) ± ly 

The pressure on the wall (for one degree of freedom) in the fast wall limit is therefore 

AP,„ = J<l>dz^J ^^^-^ \,MiE-.P.)^iy (27) 

2.7 Solutions to the hydrodynamical equations 

Detonations and deflagrations 

As mentioned, hydrodynamical treatments mostly model the plasma as a perfect 
ultrarelativistic fluid and solve the corresponding stress-energy conservation equa- 
tions. When the bubble becomes large enough and the bubble wall propagates at a 
constant velocity in relation to the medium, the bubble profile becomes a similarity 
solution to the equations, expanding linearly with time but maintaining its relative 
shape. The medium is assumed to be at rest both far ahead of and far behind the 
advancing bubble wall. This assumption dictates that no single, planar front may 
satisfy the stress-energy conservation equations (see e.g. [S])- Two possibilities are 
usually assumed to exist depending on whether the bubble wall advances at a ve- 
locity larger or smaller than the speed of sound in the medium c^: The first case 
is generally known as a detonation, the second as a deflagration. In a detonation 
no information propagates ahead of the moving wall, which hits the medium while 
it is still at rest. The medium is heated up and accelerated by the bubble wall, 
then cools down and comes to rest in a rarefaction wave while already in the broken 
symmetry phase behind the wall. In a deflagration, a shock front propagates at or 
close to the speed of sound ahead of the bubble wall. It heats up and accelerates 
the fluid in the symmetric phase before the phase transition front proper brings it 
back to rest. The common assumption of modelling the phase transition front and 
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the shock front (where it exists) as planar surfaces leads to constant values of the 
dynamic variables (fluid velocity and temperature) in the regions away from the 
fronts, including that between the phase transition front and the shock front in a 
deflagration. We will refine our calculation (for a deflagration) by splitting the bub- 
ble profile into three regions and writing the conservation equations separately for 
each. We adopt the planar approach to integrate the conservation equations across 
the bubble wall that separates the symmetric and broken phases, and across the 
shock front that propagates ahead of the wall. However (as in e.g. [IB]) we will take 
into account the sphericity of the bubble to calculate the evolution of the dynamic 
variables in the region between the fronts. In this way we may produce a whole 
profile of the bubble. We match our results to what we take to be the temperature 
of the universe outside the bubble (ahead of the shock front for a deflagration, or 
the bubble wall for a detonation). We carry out the hydrodynamical calculations 
for steady-state solutions. Finally, we investigate the possibility of runaway walls. 

Hybrid solutions 

A more careful analysis (see e.g. j?6]) reveals the possibility of so-called hybrid 
steady-state solutions in between detonations and deflagrations. Hybrid solutions 
have both a shock front propagating ahead of the phase transition front proper and 
a rarefaction wave following the phase transition front. For a given physical situa- 
tion they appear naturally with increasing steady-state wall velocity (arising with 
decreasing friction). A subsonic wall, as mentioned, propagates as a deflagration. If 
the wall becomes supersonic, however, it does not instantly become a detonation. 
Instead, a rarefaction wave appears behind the wall while a supersonic shock front 
still propagates ahead of the phase transition front. If the steady-state wall velocity 
becomes even larger, it eventually catches up with the shock front, which disappears, 
the solution becoming a pure detonation. In this work we do not concern ourselves 
with the dynamics of a rarefaction wave and therefore do not study in any detail the 
possibility of hybrid solutions. However, as noted below, our results are consistent 
with the appearance of hybrid solutions within the parameter space or our model. 

2.8 Dependence of the pressure on the wall velocity 

It has recently been established that the friction, which slows down the wall due to 
interactions between the wall and the particles in the plasma (and which exactly 
balances out with the pressure difference for a steady-state wall, see eq. ffTTjl ) reaches 
a finite limit as — > 1 [39j. This opens up the possibility of runaway walls if the 
pressure difference is higher than this limit to the friction. In the ultrarelativistic 
limit both our 'slow wall' and 'fast wall' limits to the friction integral (equal to 
the steady-state pressure difference) converge to the same integral dependence as 
E which tends indeed to a finite limit with approaching 1. It must be noted, 
however, that the 'slow wall' limit of the friction integral contains the relaxation time 
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r which should not remain constant but grow as — > 1. Assuming a dependence 
r = ro7 with tq constant the friction in this hmit diverges as 7. No such divergence 
occurs in the 'fast wall' limit. We assume therefore that our description is adequate 
in each limit. We will use this fact to link the relaxation time approximation with 
the hydrodynamical calculation of the wall velocity and produce a usable model for 
the friction in section 4. 



3 Hydrodynamical treatment 
3.1 The conservation equations 

We essentially follow the same strategy as in |12], originally adapted from [35]. We 
derive our hydrodynamical equations from the conservation of the energy-momentum 
tensor. 

The energy-momentum tensor of the system consists of that of the electroweak 
plasma (modelled as a perfect relativistic fluid) plus that of the background Higgs 
field, 



^ ~ -'■ field "T fluid ~ 



2 

+ {p + P)u^u'' - Pgf"". (28) 

As the equation of state for the system we adopt P = Pr — V{(f), T) with the 
radiation pressure P^. = aT^ = ^g*T'^, g* being the number of effective degrees of 
freedom at the temperature T. We employ the usual thermodynamic relations 



P = -f (29) 

P = f-T%^ (30) 

/ denoting the free energy per unit volume. With these the total energy-momentum 
tensor may be expressed as 



+ (^Ur - ^ ^^^'^^ ) ^'^"^ - a"" (Pr - Vi<P, T)) (31) 

where the radiative enthalpy is Ur = 4aT^. Even though this is conserved as a whole 
[dpT^'^ = 0), the fluid and field components are not: 
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We now split the conservation equations into two parts (as in [35j), making each 
part equal to plus or minus an appropriately chosen friction term which we choose 
to make dependent on a dimensionless friction parameter. Our original choice (later 
reconsidered, see I3.2p was inspired by the form of the equation of motion for the 
Higgs VEV in the 'slow wall' limit of the relaxation time approximation (12T]) . given 
the fermionic and bosonic mass dependence on the Higgs VEV. The conservation 
equations (in covariant form) may be split as 



where Vo(0) is the part of the Higgs effective potential dependent only on the Higgs 
VEV while Vi(0, T) is the part dependent also on temperature, rj is our friction 
parameter characterising the resistance of the plasma to the wall's movement and 
Tgi is a reference temperature, which we choose to be that of the plasma in the 
symmetric phase just ahead of the wall. To study the dynamics of the variables 
across the bubble and shock fronts we restrict ourselves to the 1 + 1 dimensional 
case (the spatial dimension taken as perpendicular to the front, along the direction 
of bubble expansion). We assume steady-state solutions characterised by constant 
bubble wall and shock front (where present) velocities along the spatial direction. 
Following [33] closely we arrive at the system (in the bubble wall's rest frame) 



dx^ d(j) Tsi dx 

(iar - T^-^hh. = Q (35) 

(4„T' - T'^^^^hV + P. - V(^,T) + = C, (36) 

where v and T in addition to are functions of the spatial coordinate and 7 is the 
usual relativistic factor (1 — f ^)~^/^. The value of the two arbitrary constants Ci 
and C2 is easiest to calculate in the symmetric phase, where V as well as and its 
spatial derivatives vanish. We integrate equation with boundary conditions 



0'(x) = (at both ends of the integration interval) (37) 
= (in the symmetric phase) (38) 
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at the same time solving for the velocity and temperature profiles across the bubble 
wall. We start from values of the velocity and temperature of the cosmic fluid in 
the symmetric phase ahead of the wall, Vsi and Tgi, and reach values Vh, Tb in 
the broken symmetry phase, behind the wall, see Fig [31 These values are always 
expressed in the wall rest frame. These calculations can be done by a variety of 
methods e.g. through a linearisation procedure described below. 



The shock front 

The shock front (for a deflagration) propagates across the symmetric phase so both 
the Higgs VEV and the effective potential vanish. In the rest frame of the front the 
stress-energy conservation equations simplify to 



9,[(P + p)7\] = 
9.[(P + p)7V + P] = 



which integrate to 




(39) 
(40) 

The subscripts u, s2 refer to values in the symmetric phase respectively ahead of 
and behind the shock front. 



The frame of the fluid 

For a subsonic, pure deflagration bubbl^, the velocities found in the rest frames of 
each front are related to the velocities in the rest frame of the fluid far away from 
the wall via the corresponding relativistic transformations: 



^shock front 
^^bubblc wall = 1'b 

^^fluid = (41) 

1 - Vy,Vsl 

■^^fiuid being the velocity of the fluid in the 'frame of the universe' (in which the fluid 
far ahead of and behind the bubble wall and shock front is at rest) just ahead of the 
bubble wall. 



^Note once again that we are not considering hybrid solutions. 
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Figure 3: Higgs VEV 0, velocity and temperature profiles across the (defiagration) 
bubble wall for the 0^ model with M = 800 GeV, mh = 115 GeV at the temperature 
of the universe (ahead of the shock front) = 105.49 GeV and a value of the friction 
coefficient rj = 0.398. The broken symmetry phase is on the right and the bubble 
wall propagates from right to left. The profile has been obtained solving the system 
( 1M|) -( 136|) via a linearisation procedure. 
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For a supersonic, pure detonation bubble which hits the medium while at rest 
(with no shock front) it is enough to integrate the system across the bubble 

wall, which gives us the wall velocity straight away: 



■ybubblewall — ^^sl- (42) 

Intermediate region in a deflagration bubble 

Taking the sphericity of a subsonic bubble into account requires integrating the 
conservation equations across the region between the bubble wall and the shock 
front. As we have seen, from the integration across the bubble wall we know the 
wall velocity. We now transform to the 'frame of the universe' and take advantage 
of the similarity properties of the expanding bubble to make all variables dependent 
on the single coordinate ^ = r/t. The stress-energy equations become [16] 

^~'^'-2^ + ^(l-f(e--)) 



I —v^ dP \ ^"^ 

7^ = 7 (^ - v)^- 

w at, 

Isolating derivatives, one obtains 

~ 3e(i-<)^i-< ~ m-v) 

dv 2v{l—v'^)^ — v 



d^ 3e(e-^)^i-< s{^-v) 



(43) 
(44) 



We integrate along the radial direction, starting at the bubble wall and making 
an initial guess for ^ at the shock front. Our guess is the shock front velocity so 
we can transform the final value of the fluid velocity Vs2 (just behind the front) to 
the frame of the front and use fHOj) to find the fluid velocity ahead of the front. 
We take guesses until this matches ^ for the front. We may finally calculate the 
temperature of the undisturbed universe via (!39l) . It is worth mentioning that 
for a weak phase transition the 'leap' of the variables across the shock front becomes 
extremely small (SH]- In addition, because the calculation in the intermediate region 
is done in the 'frame of the universe' where the fluid velocity is very small (unless 
the phase transition is extremely strong) we may linearise and integrate fH3|l and 
(jH]) for small v, arriving at 

_2f(eo)eo .J l^ 

T(0 =T(^o)e ^' (45) 
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where = "^bubbiewaii and f(^o) = "f^flmd- In this manner, provided we know 77, it is 
possible to look for a self-consistent, steady-state solution to the hydrodynamical 
equations for which the temperature outside the bubble is whichever T„ we wish to 
impose (the simplest assumption being to equate with the nucleation temperature 



3.2 The friction parameter and the wall velocity 

In the above analysis 77 is a Lorentz scalar which we take to be a number parametris- 
ing friction. This is what the covariant construction of the energy-momentum con- 
servation equations suggests. However this leads to immediate trouble. Note that 
the same reasoning that took us to (ITTi) . if applied to (^^, leads to a pressure dif- 
ference divergent as — > 1 because of the 7 factor present. This may not be very 
important if the bubble wall is very subsonic (as it is in supersymmetric light stop 
scenario [221 113 HH SS]) but for fast walls (or in order to investigate the possibility 
of runaway) this difficulty cannot be ignored. We are therefore forced to consider a 
more complex relationship between the friction parameter and the wall velocity. The 
simplest assumption that we can make is that the friction parameter depends on 
as ri{vuj) = Vo1~^, with rjo constant and 7 = (1 — f^)"2 ^ so that the 7 factor in (134]) 
is effectively not there. This is the assumption we shall make in our calculations. 

4 Modelling the friction 

4.1 Friction in the relaxation time approximation 

We wish to produce a reliable estimate for the friction parameter in flM]) in order to 
solve for the bubble profile and calculate the wall velocity for our model. The form 
of the friction parameter supplied by the relaxation time approximation in its 'slow 
wall' limit can be written in analogy to f l5^ through the relevant mass couplings. 



for one degree of freedom. In the last step we assumed that m is proportional to 
0. The relevant prefactors coming from r and the mass couplings depend on the 
model considered. Rather than carrying out a first-principles calculation, our goal 
is to develop a technique to relate the existing results from [SB] (arrived at using the 
full form of the collision integral) to any Standard Model-like situation (in which 





(47) 
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Figure 4: Velocity and temperature profiles across the intermediate region between 
the bubble wall and the shock front for the same deflagration bubble as in figure |3] 
with M = 800 GeV, ruh = 115 GeV at T„ = 105.49 GeV. Here the position of the 
bubble wall is on the left end of the integration interval and that of the shock front 
on the right. If the sphericity of the bubble is neglected (and a less realistic planar 
approximation adopted instead), v and T do not vary across this region. 
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friction is carried by the top and weak gauge boson degrees of freedom). A similar 
analysis for the case of friction dominated by light stops was carried out in 

Note that we are free to change integration variables in fH7|) from p ^ p' = ^. 
Since E = a/ m'^ + the friction term suggested by the relaxation time approxima- 
tion in this form becomes 

dm^ f d^p 

-df{p,x) ~ 



(27r)32E 



(J)'^ (f)' T'JV 



e7(V(f)'+p"-^) 

(27r)3 ((^)2 +p'2) ^^^(^(^jF^^^iiZ.) ^ ^^2' 

(4J 



Given the general dependence of the mass the integral contains factors of ^. This 
strongly suggests that friction is dependent on the strength of the phase transition, 
commonly expressed in the literature by the parameter = f^, that is, ^ evaluated 
at the critical temperature. If we wish to approximate the friction integral by a 
single numerical contribution (from each degree of freedom), one way would be to 
approximate the ^ factors by ^ (or an appropriate related value like | |f, the 

value of ^ roughly half way across the wall at Tc) and integrate over momentunj^. 
Alternatively, one can plot the spatial dependence of the friction term itself (peaked 
because of the (f)' factor, which goes to zero away from the bubble wall on either side) 
and replace the integral by a constant which equalises the peaks (Figure E]). This can 
be done, for example, through the 2-parameter hyperbolic tangent Ansatz commonly 
used to approximate the shape of the wall (see e.g. p^). The variation of the Higgs 
VEV across the bubble wall (assumed planar) is given in this approximation by (Fig 
ED 

0(z) = y (l + tanh |) (49) 

where L is the wall thickness and 0o the value of the Higgs VEV in the broken 
symmetry phase behind the bubble wall as read from the effective potential. As the 
Ansatz is written in the rest frame of the advancing steady-state bubble wall there 
is no time dependence. A quick way to obtain a usable value for L is through 

- I (50) 

with Vfc the height of the potential barrier between the two minima of the effective 
potential (this is suggested by the simplified 2-parameter scalar potential ^(0) = 
|0^(0 — 0o)^ with the equation of motion ^ = into which the hyperbolic 

tangent Ansatz may be substituted). 



®Note that there are other choices as regards expressing the strength of the phase transition, 
namely as a function of the nucleation (^„ = |?^) or broken symmetry phase (^6 = 1^) temperatures. 
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Figure 5: Spatially varying part of the friction term according to the relaxation 
time approximation (slow wall hmit), 0V / 12-j,y{2E)^ (^Jj(e-v^Z7±i)2 (solid line), and 
the same term with the momentum integral replaced by a fitted constant (dotted 
line), C 4>'^4'' , for fermions (left) and bosons (right) in an example case. We use the 
hyperbolic tangent Ansatz to approximate the bubble profile and assume 0o = 100 
GeV in the broken symmetry phase, T = 100 GeV, Ly^ ■ T = 15, and a mass 
dependence m = -^(p- 

Using either of these approaches, or a comparable simplification, the relaxation 
time approximation model supplies a numerical coefficient for a description of the 
friction of the form rjivw) = Vol~^- Because we assume that rjo has no dependence 
on Vw we set = in the momentum integral in f HHj) . 

4.2 Fixing the friction parameter 

As mentioned, the prefactors for the friction integral for each massive degree of free- 
dom in the plasma that couples strongly to the Higgs can be calibrated from existing 
results, in our case those for the Standard Model in |36]. This is because friction 
processes are a function of the particle content of the model and the relevant particle 
interactions, and therefore, as long as we remain within a Standard Model-like sit- 
uation, the prefactors to the momentum integrals (the relevant coupling constants 
and the chosen form for the relaxation time approximation) should hold. This is, 
for example, the case in the dimension-6 extension of the SM studied below. 

The standard model 2-loop, high-temperature thermal effective potential used 
in [SB] is 

Feff (0, T) = D{T' - T^)<P' - CT'^' ln(|) - ET^' + (51) 
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Figure 6: Sample hyperbolic tangent Ansatz for the Higgs VEV across the bubble 
wall. The broken symmetry phase is on the right hand side with an assumed VEV 
of 00 = 100 GeV. We take L = 50 GeY'\ 
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with 



1 
1 



8v^ 



C = ^{lA2gt + 4.8glX-6X') 
E = -L[Ai^Y + 2i^f + i3 + 3'-')X''] 

l/TT Vq Vq 

B = ^(2mt + m^-4m^) 



^0 = \l^iml-8Bvi) 



and 



2mw 



= 80.4 GeV 
m, = 91.2 GeV 
mt = 174.0 GeV 
vo = 246.0 GeV 
ttb = 49.78019250 
ttf = 3.111262032. 

(52) 

We show the values of rjo which reproduce the values of Vw in [36] in table [1] 
(denoted by tiq^sm)- We now wish, as mentioned, to produce a prediction for the 
friction coefficient for our model. The authors of estimate 60% of the Standard 
Model friction to come from fermions and 40% from bosons. We call the values of the 
bosonic and fermionic numerical coefficients obtained from the momentum integrals 
in eq f HHjl for each case in [36] (and characterised by a value of the strength of the 
phase transition ^ = ^) lobiO^ ^ofiO- With this calibration, the friction parameter 
for SM-like friction, e.g. within the dimension-6 model (characterised by a specific 
value of the strength of the phase transition ^) may be found as 

^o(0 = Vo,SM (0.6/41^ + 0-4/#y) • (53) 



5 The dimension-6 extension to the SM 

It has been suggested [50], EH [52] that new physics may appear as higher- dimensional, 
nonrenormalisable operators added to the Standard Model Higgs scalar potential. 
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Table 1: Values of the Standard Model quartic coupling A, Higgs mass, strength of 
the phase transition f„ = nucleation temperature and fitted friction coefficient 
7] for relevant cases in [36] . 



Ay 


m;,(GeV) 


f _ t'n 


t;^^, (from [21 


3) T„(GeV) 


^70,5^/ (fitter 


0.023 





0.98 


0.374 


57.192 


0.5628 


0.030 


50 


0.80 


0.392 


83.426 


0.5286 


0.041 


68 


0.65 


0.412 


100.352 


0.6431 


0.050 


79 


0.58 


0.428 


111.480 


0.6705 


0.060 


88 


0.53 


0.441 


120.934 


0.6707 



getting around the constraints posed by present experimental bounds on the Higgs 
mass and providing additional sources of CP violation. The addition of dimension-6 
operators to the Higgs potential |33l |32l [53] has been proposed, numerical calcu- 
lations having shown that further higher-order terms suppressed by the same low 
cut-off scale as the dimension-6 terms give corrections of only a few percent to the 
strength of the phase transition |f [32]. The dynamics of the electroweak phase 
transition in this scenario are parametrised by the Higgs boson mass ttih and the 
cut-off scale M. In this situation the quartic coupling of the Higgs potential may as- 
sume negative values. Such dimension-6 operators could stem from integrating out 
a massive degree of freedom like a scalar singlet [32], or alternatively from sources 
such as strongly coupled gravity [33] . New physics with a comparatively low cut-off 
scale may lead to non-standard signals which could be detected in the near future, 
such as modified Higgs self-couplings [S2] . 



5.1 The effective potential 

As dimension-6 effective potential we take 



8M2 



+ 20^T' + 0'T4) (54) 



where Q = mtop and cp ~ 13.94. M and rrih are the free parameters of the model, /i 
and A are found through the conditions for the zero-temperature potential {vq = 246 
GeV) 

^,.,,.0, ^ . A,. , _ JL,.,.,„ (1^) ,55) 
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which take the form 



dVes{(p,0) 
(90 



<p=vo 







502 



(p=V0 




(56) 



6 Results for the dimension-6 extension 

6.1 General solutions to the hydrodynamic equations. 

We begin by studying the dependence of the steady-state wall velocity for pure 
deflagration and detonations solutions (excluding hybrids) on arbitrary values of 
the friction parameter for a specific choice of dimension-6 model parameters. We 
present the results of such an analysis for an arbitrary choice of model parameters in 
Figure [71 expressing the wall velocity as a function of the temperature of the universe 
Tu for values of ?7 = 0.3, 0.4, 0.5. The general shape of our solutions agrees with 
previous studies (see e.g. |1Q]). At high values of the temperature of the universe 
only subsonic solutions are allowed. The wall velocity increases with decreasing Tu 
and is obviously higher for lower values of t], becoming larger if we choose to neglect 
the sphericity of the expanding bubble and assume planar symmetry in the region 
between the shock front and the bubble walfl. As Tu keeps decreasing two branches 
of additional, supersonic solutions appear. Of these only the upper branch is physical 
(note that the shape of the lower branches implies that Vw would decrease with lower 
?7, and again decrease as T„ decreases and the free energy difference released by the 
phase transition becomes larger). 

A related issue is the stability of the expanding bubbles, examined in ^7]. That 
reference concludes that subsonic bubbles may become unstable to the appearance 
of perturbations larger than a critical size provided that, essentially, the wall velocity 
decreases with decreasing Tgi. We may equivalently plot T„ vs T,! (Figure [8]), finding 
the values of T„ for each value of the friction coefficient below which subsonic bubbles 
may become unstable. For the choice of model parameters in Figure [TJ and as a 
general rule, this is roughly the point at which supersonic solutions first appear, 
and also the region in which hybrid solutions would naturally exist. Low- velocity 
subsonic solutions for high Tu are always stable. 

6.2 The wall velocity in realistic cases. 

We conclude by producing predictions for the numerical value of the friction parame- 
ter in a realistic setting and, through our full hydrodynamic calculation, for the wall 
velocity. Relying on recent LHC results we do this for a fixed Higgs mass rriH = 125 
GeV and study the variation of the friction and v^, with the cut-off scale M. We 
wish to investigate the impact on our studies of the criterion for runaway walls from 

^As noted, with such a simphfication the dynamic variables v, T do not change across this 
region. 



25 




Figure 7: Steady-state bubble wall velocity vs temperature of the universe for three 
values of the friction coefficient 77, 0.3 (blue), 0,4 (green) and 0.5 (yellow) for the 
dimension-6 model with M = 800 GeV, rrih = 120 GeV. The continuous lines below 
the horizontal line which marks the speed of sound are subsonic solutions, the dotted 
lines marking solutions found through neglecting the sphericity of the bubble. At the 
horizontal error bars the stability criterium for subsonic solutions changes sign (the 
stability region lies to the right of the mark). The crosses above the horizontal lines 
indicate supersonic solutions in which the bubble wall hits the medium at rest. The 
two vertical lines mark the nucleation and finalisation temperatures for the phase 
transition for these parameters. Note the two branches of supersonic solutions for 
each 7]. In this example and for these values of rj supersonic steady-state solutions 
would be excluded but stable subsonic ones allowed throughout the duration of the 
phase transition. 
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Figure 8: Temperatures in the symmetric phase just ahead of the bubble wall vs 
plasma temperature for the parameters of figure ([7j) and values of the friction pa- 
rameter rj = 0.3,0.4,0.5 (blue, green and yellow). Right of the minimum T^i for 
each r] is the stabihty region for subsonic solutions according to the criterion in [37] . 
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Table 2: Values of the friction coefficient and the wall velocity for the dimension-6 
extension to the Standard Model, as a function of the cut-off scale M. Two sets of 
values of the friction coefficient are given, both calculated through the relaxation 
time approximation based on: 1) Existing Standard Model wall velocity calcula- 
tions (?7sm); and 2) The runaway criterion in ^39j [f]^). The wall velocity is calculated 
through the Standard Model-based value of the friction parameter, r/sM- An approx- 
imate value of Vw for r^^unaway is also shown. We assume a Higgs mass mh = 125 GeV. 
For each value of M the nucleation temperature and the corresponding strength of 
the phase transition S,n = p- are given. 



M 


T 


e = 0o/T„ 


'7SM 








900 


130.08 


0.87 


0.615 


0.674 


0.312 


0.284 


800 


120.57 


1.30 


0.458 


0.487 


0.366 


0.344 


700 


106.28 


1.86 


0.340 


0.352 


0.447 


0.432 


600 


80.06 


2.90 


0.203 


0.208 


0.829 


0.809 


590 


75.59 


3.11 


0.184 


0.189 


0.920 


0.898 



[39], which establishes that, if the slowing effect of the plasma on the advancing wall 
in the limit f^^, — )■ 1 is not sufficient to counter the accelerating force of the vacuum 
potential, the wall will run away. The phase transition becomes stronger, and the 
bubble wall faster, with decreasing M. Correspondingly we find that (for a Higgs 
mass of 125 GeV) the runaway criterion is satisfied for M < 582 GeV and a strength 
of the phase transition in = ^ 3.3. 

The runaway criterion also provides us with an alternative calibration point for 
the friction parameter. Instead of relying on values of the parameter fitted to existing 
Standard Model wall velocity calculations (as we wrote in eq. f l53p ) we may calculate 
the value of the friction parameter T^runaway that gives v& Vw ^ I for M ~ 582 GeV, 
niH = 125 GeV, and write the calibration formula based on that /^runaway and the 
corresponding value of the strength of the phase transition. In Table [2] we present 
the values of the friction parameter calculated on the basis of both calibrations 
(from existing Standard Model results and from our runaway point), which show 
remarkable consistency. We provide the values of the wall velocity calculated on 
the basis of our original calibration (based on the results in [36]) for a range of 
M values. The same results for the wall velocity are plotted in figure |9] alongside, 
for comparison, approximate values for found using the corresponding friction 
parameter values from the runaway point calibration (These approximate velocities 
are found by taking advantage of the fact that calibrated friction parameter values 
and wall velocities show a rough correspondence const ant ) . Figure [10] shows 

the calculated wall velocities as a function of the strength of the phase transition at 
the nucleation temperature ^n- We want to stress that our linearized treatment of 
the hydrodynamic equations breaks down close to the speed of sound. As a result 
no cases in Figure dU] have wall velocities close to the speed of sound. 
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Figure 9: Wall velocity in the dimension-6 model as a function of the cutoff scale 
M for mn = 125 GeV calculated through the friction parameter obtained from the 
Standard Model calibration (in black, above) and the alternative calibration from 
the runaway criterion in [39] (in blue, below). The horizontal line marks the speed 
of sound in the medium. Wall velocities predicted in this way become supersonic 
below approximately M ^ 640 GeV. Note, however, that our linearized treatment 
of the hydrodynamic equations breaks down close to the speed of sound. 
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Figure 10: Wall velocity in the dimension-6 model for mn = 125 GeV as a function 
of the strength of the phase transition at the nucleation temperature ^„ = |r-. 



30 



7 Conclusions 



The bubble expansion velocity is a crucial ingredient in any computation of relic 
signals from a first-order electroweak phase transition. For instance, standard elec- 
troweak baryogenesis relies on a subsonic wall to allow for efficient diffusion of chiral 
charges into the symmetric phas^. Strong gravitational waves signals, on the other 
hand, require a supersonic bubble wall. 

The microscopic determination of the wall velocity, even in the semiclassical 
approximation, is a very difficult task, which requires the simultaneous solution 
of the Higgs equation of motion and Boltzmann equations describing the out of 
equilibrium plasma. To date this has been achieved only for the cases of the SM 
[55] and the light stop scenario of the MSSM [55] . 

The value of the approach developed in this work is that it uses the results of 
these microscopic computations and allows to transfer them to related scenarios, 
where the Higgs dynamics is different but friction remains essentially the same as in 
a) [2S] (dominated by tops and weak gauge bosons) or b) [3H] (dominated by light 
stops). Other prominent extensions of the SM, such as the Two-Higgs-doublet model 
or SUSY models without light stops should, for instance, be quite well described by 
case a). 

The computation of the wall velocity, as laid out in this work, consists of series 
of steps: Firstly, the temperature of the phase transition has to be determined as 
described in section 2.2. The main equations to consistently solve then are (34) - 
(36), which describe the Higgs field coupled to an ideal fluid characterised by a fluid 
velocity and temperature. The unknown friction parameter rj can be obtained from 
eq. (53), where the momentum integrals If,biO devive from eq. (48). For the SM 
one has r^o.SM ~ 0.6 (see Table 1). In the case of a deflagration, eqs. (45,46) have to 
be used to connect the bubble wall solution to the shock front. 

In the last part of this work we have applied our approach to the SM augmented 
with a 4>^ operator in the Higgs potential. For ^ ~ 1, i.e. a phase transition strong 
enough to prevent baryon number washout, we flnd a wall velocity = 0.3 — 0.4, 
i.e. subsonic walls, which support the standard picture of electroweak baryogenesis 
assumed in [33]. As the strength of the dimension-6 operator increases, the phase 
transition gets stronger and the walls become faster and at some point become 
supersonic (see Fig. 9). Close to the speed of sound our computation is not reliable, 
as in the evaluations we linearize in the plasma velocity. For very strong phase 
transitions, ^ > 3.3 we observe runaway behaviour of the bubble wall, consistent 
with [M]. 

It is interesting to note that ref. [30] provides a direct criterion for runway be- 
haviour. We can use this to flx rjo^sM independently of [36]. The two determinations 
of the friction parameter agree well, which is very satisfactory (see Table 2). We 
take this as additional support for the method presented in this work. 

It will be very interesting and fruitful to use this formalism to compute the wall 

®See, however, ref. [54] for an interesting alternative, 
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velocities in other promising extensions of the SM, such as the Two-Higgs-doublet 
model or non-minimal SUSY models to arrive at more reliable predictions for relics 
of a first-order electroweak phase transition. 
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